A general test of the Copernican Principle 
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The recent discovery of apparent cosmic acceleration has highlighted the depth of our ignorance of 
the fundamental properties of nature. It is commonly assumed that the explanation for acceleration 
must come from a new form of energy dominating the cosmos - dark energy - or a modification 
of Einstein's theory of Relativity. It is often overlooked, however, that a currently viable alterna- 
tive explanation of the data is radial inhomogeneity which alters the Hubble diagram without any 
acceleration. This explanation is often ignored for two reasons: radial inhomogeneity significantly 
complicates analysis and predictions, and so the full details have not been investigated; and it is a 
philosophically highly controversial idea, revoking as it does the long-held Copernican Principle. To 
date, there has not been a general way of determining the validity if the Copernican Principle - that 
we live at a typical position in the universe - significantly weakening the foundations of cosmology 
as a scientific endeavor [J]. Here we present an observational test for the Copernican assumption 
which can be automatically implemented while we search for dark energy in the coming decade. Our 
test is entirely independent of any model for dark energy or theory of gravity and thereby represents 
a model-independent test of the Copernican Principle. 
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Introduction The Copernican Principle is at the root 
of cosmology, having played a pivotal role from the be- 
ginning of the modern history of the subject. Einstein 
erroneously imposed the Copernican Principle in both 
space and time by demanding a static universe, forcing 
him to introduce the repulsive cosmological constant. It 
is ironic that this very same cosmological constant may 
be responsible for the apparent accelerated expansion of 
the universe, thereby mysteriously placing the beginning 
of life on earth disconcertingly near the special cosmic 
time when this acceleration began. The magnitude of the 
cosmological constant required to achieve this is inexpli- 
cably small given our current understanding of quantum 
physics. Given these problems, might it not be just as 
natural for us to be living at a special place in a universe 
with no cosmological constant or dark energy, despite the 
charged religious connotations it would imply? Although 
violation of the Copernican Principle would not solve the 
cosmological constant problem, it would allow explana- 
tion of the data without the need for any dark energy. Of 
course, the bedrock of the standard model relies on the 
reasonable assumption that we are not living at a special 
place. But to place cosmology on a firm scientific ba- 
sis it is ever more clear that this matter must be settled 
observationally without theoretical bias. 

Until now this meant trying to develop inhomogeneous 
and/or anisotropic models to the level of sophistication 
that they could be tested against observations - an 
extremely laborious task given that the Einstein Field 
Equations are nonlinear coupled partial differential equa- 
tions. The usual argument for ignoring these difficult 
models comes from Occam's razor: why bother since we 
have a simpler model which is working so well? While 
radial inhomogeneity can fit any Hubble diagram per- 
fectly 0, 0, [1| , in general it requires the introduction of 



an infinite number of free parameters to do so, clearly 
disfavored from a Bayesian point of view. However, it is 
not clear that one- or two-parameter models of radial in- 
homogeneity would not give as good a fit as homogeneous 
dark energy models with the same number of degrees of 
freedom, which leaves us in the same position of arbitrary 
model building that we have at present in the standard 
paradigm which instead parametrises an otherwise arbi- 
trary equation of state for the dark energy, w(z). 

Instead what we need is a model-independent test to 
determine if the Copernican assumption is in fact vio- 
lated. One elegant suggestion to falsify the CP uses ob- 
servations of the CMB from inside our light cone 0,[3|, ei- 
ther by testing for deviations from a blackbody spectrum, 
or deviations from isotropy at distant clusters. These 
important methods are very difficult and require precise 
control of systematic effects. Furthermore, these tests 
have a fundamental limitation. Even if we were to deter- 
mine isotropy of the CMB for all observers in the uni- 
verse it would still not be enough to prove homogeneity, 
unless we also proved that all fluid components of the 
cosmic energy budget are co-moving, barotropic perfect 
fluids 0. 

We present here a new test for the CP which relies 
on a consistency relation that exists within the homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models between luminosity or area dis- 
tance and the Hubble rate, both as a function of redshift. 
We show that this is not satisfied for radially inhomoge- 
neous models in general, providing just such a mechanism 
to test whether we live in an approximation to an FLRW 
universe or, instead, near the center of a spherically sym- 
metric one. 
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Consistency Relation In FLRW models, the luminos- 
ity distance may be written as (in units where c = 1) 
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where fifc is the curvature parameter today, and the ex- 
pansion rate H(z) is given by the Friedmann equation 
and takes the value H = H(0) today. The area distance 
is defined using di, = (1 + z) 2 d^ and another distance 
measure we will find useful is D = (1 + z)dA- We may 
rearrange Eq. ([I]) to give an expression for the curvature 
parameter in terms of H (z) and D(z): 



[H{z)D'{z)] 2 - 
[H D{z)f 



(2) 



where ' = d/dz. This tells us how to measure the curva- 
ture, in principle, from distance and Hubble rate obser- 
vations, independently of any other model parameters or 
dark energy model. Remarkably this tells us the curva- 
ture today from these measurements at any single red- 
shift. 

Given that the curvature parameter is independent of 
redshift, we may differentiate this to obtain an expression 
which must equal zero. The factor responsible for this is 



(z) = 1 + H 2 (DD" - D' 2 ) + HH'DD' 



(3) 



which must be zero in any FLRW model at all redshifts, 
by virtue of Eq. (TT]) . This function may also be derived by 
equating the two reconstructed functions for w(z) given 
by Eqs. (2) and (3) in [|. 

As we have not utilised the Friedmann equation, the 
derivation of ^(z) relies only on the metric of spacctime, 
and not on the theory of gravity, nor on any matter model 
present - it is therefore a purely geometric function. Con- 
sequently, if the FLRW models are indeed the correct 
background model, then we should expect to measure 
c &{z) w (up to the amplitude of perturbations) in the 
real universe at all redshifts. 

We may estimate the errors on ^(z) using series ex- 
pansions for D(z) and H(z), but with different FLRW 
parameter values, and dark energy equations of state. 
At leading order we have 
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i.e., the difference between the separate deceleration pa- 
rameters measured using D and H. Hence, at low red- 
shift, our error on 'jf(z) is roughly the same as go a s 
estimated using H(z). 

A critical aspect of our test is that H(z) measure- 
ments must come from a different data set than used 
for distance data. Presently most estimates of H(z) 
rely implicitly on Eq. ^ or some other aspect of the 
FLRW models themselves. However, passively evolving 



objects such as luminous red galaxies provide a family 
of roughly synchronised clocks at moderate to high red- 
shift, from which we can measure H{z) directly using 
H(z) = —1/(1 + z)i'(z) [IH. Alternatively, measuring 
the change of redshift, z(z), of an object over a period 
of several years can also provide H(z) [nj, and so be 
used as a probe of the Copernican Principle based on the 
arguments presented here [lij ]. 

Prospects for constraining c <o(z) in the near future 
as an independent confirmation of the standard model 
are exciting. For example, finding w(z) from distance 
data requires knowledge of D"(z), while reconstructing 
it from H(z) alone requires knowing H'(z) Tak- 
ing derivatives of observable functions introduces extra 
errors, so it might in practice be simpler to measure 
{[H(z)D'(z)f - l}/[H Q D(z)} 2 at different redshifts, and 
check that it yields the same value as it should by virtue 
of Eq. ([2]). In FLRW models this tells us the value of 
the curvature today; in more general models this quan- 
tity will not be as simple as this. This is an important 
consistency check on the FLRW models. If it is found 
to change with redshift then that is equivalent to finding 
c (o{z) ^ [l[. Thus, while we try to search for w(z), we 
can measure ^(z) at the same time to a similar degree 
of accuracy. 

But what would a measurement ^ imply, even 
if only at one redshift? Since it's not dark energy or de- 
viations from general relativity, the origin must be a dif- 
ferent cosmological geometry which is not homogeneous. 
Clearly, then, measuring ^(z) ^ at any redshift would 
pose serious problems for cosmology. 

We have given a necessary condition for violation of 
the Copernican Principle, but is it sufficient? In other 
words, could the principle be violated while still hav- 
ing 'rf(z) — 0? Consider Lemaitre-Tolman-Bondi (LTB) 
models, which can fit the Hubble diagram without in- 
voking dark energy 0]. These are spherically symmetric 
models with each shell at radius r evolving as a separate 
FLRW model. Clearly the generic scenario is for each 
shell to be characterised by a different and so have 
ft (z) 7^ 0. Indeed, we show in the appendix that ff(z) is 
a freely specifiable function in these models, so ^oiz) = 
which arc not FLRW is a very limited subclass within 
the LTB family. 

Alternative Tests In addition to the key test just pre- 
sented we present two other novel approaches to the prob- 
lem. Future cosmological data will allow us to measure 
w(z) (assuming FLRW and an appropriate parametri- 
sation) to high accuracy [HI]. With w(z) known, one 
can make predictions for what will be observed with the 
Alcock-Paczynski (AP) test 14{. The AP test relies on 
the idea that a spherical object in real space will appear 
oblate or prolate in redshift space due to the fact that 
the radial size, L\\ is determined by Az/H(z) while the 
transverse size, L± is determined by cf^z). By demand- 
ing that L± = Lii one can determine cosmic parameters, 
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in particular, w(z). 

The canonical AP test in modern cosmology is pro- 
vided by the Baryon Acoustic Oscillations (BAO) which 
uses the excess signal at ~ 150Mpc in the two-point cor- 
relation function as a standard ruler. By observing this 
in both the angular and radial directions BAO will pro- 
vide both cIa(z) and H(z), and hence allow us to measure 
w(z) to similar accuracy as SNIa [HI]. 

If the Copernican Principle is violated through signif- 
icant radial inhomogeneity then clustering in real space 
will not be isotropic. Hence the assumption underlying 
the AP and BAO tests (that the "bump" in the two- 
point correlation function is isotropic) will be invalid. If 
one incorrectly assumes isotropic clustering one will re- 
construct a w(z) that disagrees with that from SNIa or 
other measurements. Put another way this would ap- 
pear to violate distance duality (the equivalence between 
luminosity and angular-diameter distances) Tan- 
talisingly such a mismatch has actually been observed 
between current BAO and SNIa data at about the 2a 
level 0. 

Indeed there are at least three ways in which violation 
of the Copernican Principle would deform standard BAO 
results. First, the sound horizon - the fundamental stan- 
dard ruler of the BAO method - would be different in the 
|| and _L directions (relative to the observer) and in gen- 
eral will differ from the standard FLRW value (~ 150Mpc 
today). Both of these will cause problems. Secondly, the 
subsequent expansion of the sound horizon in the || and 
_L directions will be governed by different Hubble rates 
Hu and H± (see eq. [5]). Hence, even if the sound horizon 
were isotropic and equal to the standard FLRW value 
at decoupling, the subsequent evolution would cause bi- 
ases in measured cosmology. Finally, redshift distortions 
which need to be subtracted [l(| in order to apply the 
AP test, will be significantly more complicated to model 
and crucially will be correlated with the background ex- 
pansion, unlike the FLRW case. To see this consider a 
galaxy moving in a gravitational potential, $ = $m + 5$ 
which is the sum of the usual perturbative gravitational 
potential and the radial inhomogeneity of the LTB back- 
ground. Since Vj_$|| = 0, the || velocity of the galaxy 
will be different from the FLRW case even if the (5<E> con- 
tribution is identical in both cases. This will modify the 
Doppler shift contribution to the redshift of the galaxy 
and hence will modify the usual Kaiser and "Fingers of 
God" effects in a way that is strongly correlated with the 
precise nature of the radial inhomogeneity. 

Each of these effects will cause a bias in general in 
the estimation of dark energy parameters obtained by 
assuming a FLRW background which will manifest as 
a mismatch with the parameters derived from distance 
measurements derived from supernovae or lensing. It is 
also perhaps worth noting that the growth rates will be 
different in the || and _L directions. 

A significantly more exotic test for the Copernican 



Principle is provided by the realisation that the c2l(z) 
relation need not be single- valued; see Figs 5-8 of [4|]. 
The same redshift can correspond to more than one 
distance. Clearly in FLRW this is impossible and the 
occurance is related to the possibility of gravitational 
blueshifts/redshifts adding to the standard expansion 
redshift. Hence, object A which is further away than 
object B, may still have the same redshift if it lies at 
a suitably higher "potential energy" which compensates 
the difference in cosmic expansion. Hence, the volume as 
a function of redshift can exhibit pathological behaviour 
which would be visible in number-counts or a large dis- 
persion in SNIa distances around the same redshift pro- 
viding a clear signal of violations of the CP. 

Discussion On what timescale can we expect these 
various tests to be implemented? Future BAO surveys 
such as WFMOS will measure oIa(z) and H(z) to better 
than 2% and large radio surveys such as HSHS and SKA 
will measure them to better than 1% in redshift bins of 
width Az ~ 0.2 17, lfi| . At its simplest our test can be 



implemented by measuring Qk at two different redshifts. 
Combinations of future all-sky lensing, BAO and SNIa 
surveys will allow a measurement of with accuracy ~ 
0.04 2l| which sets the approximate scale of how well we 
will test the Copernican Principle in the next decade and 
which should allow us to strongly rule out inhomogeneity 
as the source of oversized cosmic distances. Beyond these 
surveys one can envisage large, low-frequency HI surveys 
able to measure (1a{z) and H(z) at high redshift, z ~ 
10 — 100 which will allow the testing of the Copernican 
Principle in the first billion years of cosmic history when 
most perturbation modes were still in the linear regime. 

While the Copernican assumption is unlikely to be dra- 
matically wrong in our Hubble sphere, it is possible that 
the dark energy phenomenon is actually concealing some- 
thing even more bizarre than a misconception of funda- 
mental physics. Indeed, with suggestions such as chaotic 
inflation, everlasting bouncing universes and the like, it 
is statistically possible that our part of the universe has 
some unusual features on a spatial, instead of a purely 
temporal, dimension, see e.g. 22j. This cannot be ruled 



out at the moment on anything other than philosophical 
grounds. 

We have presented a new, straightforward, test of the 
Copernican assumption which may play an important 
role in our understanding of dark energy. While we strive 
to determine the function w(z), we have shown that 
we can simultaneously constrain our Copernican func- 
tion ^(z) to a similar degree of accuracy. At any red- 
shift, a measurement of ^(z) ^ would imply that the 
FLRW models arc the wrong foundation for cosmology 
and that something more sophisticated must be consid- 
ered instead. Even if the cosmos is FLRW on average, 
measurement of ^(z) at small redshifts will allow us to 
probe the scale at which homogeneity sets in. The pro- 
posed test does not depend on any theory of gravity nor 
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on our understanding of dark energy, but relies only on 
the geometry of the FLRW models themselves. 
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lis, Pedro Ferreira, Ariel Goobar, Charles Hellaby, Roy 
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Appendix: c &{z) in LTB models The general spher- 
ically symmetric metric for an irrotational dust mat- 
ter source is the Lemaitre-Tolman-Bondi (LTB) metric, 
given in synchronous comoving coordinates as 



ds 2 



-dt 2 



[d r R{t,r)f 
1 + 2E(r) 



dr 2 + R 2 (t,r)dft 2 



where d r R(t, r) = dR{t, r)/dr, and dfl 2 = dd 2 +sin 2 
is the metric on a unit 2-sphere. The function R = R(t, r) 
is the areal radius. The function E = E(r) > —1/2 is 
an arbitrary function of the LTB model, representing the 
local curvature. 

Solving the Einstein field equations, we have a gener- 
alized Friedmann equation for the angular Hubble rate: 



H±(t,r) 



d t R(t, 
R(t,r 



= R(t,rY 



<2M(r) 
R(t,r) 



2E(r) , 



(4) 

where M(r) is another arbitrary function of the LTB 
model that gives the gravitational mass within comoving 
radius r, which may be derived given the energy den- 
sity p(t,r). There is one further free function of radius, 
the bang time. Of the three radial degrees of freedom, 
one is gauge, representing our choice of radial coordinate, 
while two are genuine physical degrees of freedom, and 
can be specified arbitrarily. We shall assume that the 
coordinate degree of freedom removes (or specifies) the 
bang time function: thus, specifying {M(r), E(r)} fully 
specifies the LTB model. 

In LTB models the expansion in the radial direction, 
H\\ (t, r), is not the same as that in the angular direction, 
but is given instead by [13] 



H\\(t,r) 



d t d r R(t, r) 
d r R{t, r) 



H±(t,r) 



d r H ± (t,r) m 
<9 r ln R{t,r) : 



(5) 



in the FLRW limit d r H±(t,r) = 0, so the two Hubble 
rates are the same. This can present different ways to 
view 'if(z) in LTB models, as we can use whichever H (z) 
we like. 

Finding the area distance explicitly in terms of redshift 
is not simple, and is discussed in detail in Q . In the usual 
procedure, one chooses a radial coordinate which flat- 
tens out the central observer's null cone. Evaluating the 
function R on the null cone then turns it into the area- 
distance as a function of redshift: (^(z) = i?| nu ii cone- 
However, it has been shown in [25| that di,(z), or equiva- 
lently dA{z) = R(z) (and so D(z)), may be considered as 



a free function of the LTB model instead of M(r). There- 
fore, the LTB model may be directly specified by the free 
functions {D(z), E(r)}. On the null cone we may find the 
radial coordinate r(z), implying that E(r(z)) = E(z) is 
a free function instead of E(r). As a free function, we 
may replace E(z) with H±(z) using Eq. ^ evaluated on 
the null cone, which means that an LTB model can in 
fact be specified by the free functions {D(z), H±(z)} or 
{D{z), H\\(z)}. 

So, as far as ^(z) is concerned, we now see that in 
an arbitrary LTB model, it can be anything one chooses, 
whether we are to interpret H(z) as Hj_(z) or Hu(z), 
or a suitable combination of them.Q If it is zero for 
both H±(z) and H\\(z), potentially giving a class of 
LTB models which would fail our Copernican test, then 
that leaves just one degree of freedom - exactly as in 
FLRW models with dark energy. If this free function is 
H(z) = H±(z) = Hu(z), say, then D{z) must be given 
by Eq (1), which can be shown by integrating ^(z) = 0. 
Although we have not shown that these models are nec- 
essarily FLRW, this would have to be a very restricted 
family within the full LTB class. It would be interesting 
to determine the exact conditions under which %?(z) = 
observed from one location is a sufficient condition for an 
expanding spacetime to be FLRW. 

This 'Copernican function' therefore must be consid- 
ered, as far as a test for the CP is concerned, as essentially 
free, and must therefore be determined by observations. 
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[1] We thank Ruth Durrer for pointing this out to us. Note 
that we must use the relation given by Eq. ((2} for this 
to work; measuring curvature at different redshifts by 
standard tests will be contaminated by our assump- 
tions about the dark energy equation of state. 

[2] We have verified that in fact ^(z) is non-zero for a 
class of LTB models, using the code presented in [23| . 



